A quantum motor: directed wavepacket motion in an optical lattice 
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We propose a method for arbitrary manipulations of a quantum wavepacket in an optical lattice 
by a suitable modulation of the lattice amplitude. A theoretical model allows to determine the 
modulation corresponding to a given wavepacket motion, so that arbitrary atomic trajectories can 
be generated. The method is immediately usable in state of the art experiments. 
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The fine manipulation of wavepackets is a fundamental 
requirement in a large number of fields. Particularly im- 
portant examples are quantum transport , "quantum 
simulators" that aim to reproduce solid state models (gI- 
12 1, quantum information 13l-[l6t and quantum metrol- 
ogy \XM' The possibility of trapping very cold atoms by 
light was a decisive step leading to enormous progress in 
this field. Tailored optical potentials, created by mul- 
tiple interfering beams interacting with cold atoms, al- 
low to trap and to guide atomic wavepackets for long 
times (compared to the dynamics of the atom external 
degrees of freedom) and distances (compared to the de 
Broglie wavelength). Moreover, using far off-resonance 
beams reduces decoherence effects to negligible levels. 
With such techniques old problems of quantum dynamics 
have been experimentally studied, as the elusive Bloch 
oscillations [H, [l^ or quantum chaos [i^-Il^- 
emerging, and rapidly developing, field of quantum in- 
formation, controlled motion in optical lattices provides 
ways to manipulate q-bits (27I, \2^. The next step in the 
development of such techniques is the shaping and dis- 
placing of wavepackets at will, and this is the problem 
addressed in the present work. We consider the motion 
of a wavepacket in a driven two-dimensional (2D) optical 
lattice (the generalization to the 3D case is straightfor- 
ward). By carefully engineering the temporal driving of 
the optical potential, we demonstrate a way to coher- 
ently impinge to an atom in an lattice almost any kind 
of motion, including coherent rotations of its wavepacket. 
Tailored wavepacket motions have been previously stud- 
ied both theoretically and experimentally by using Bloch 
oscillations jll, [2§| , but in this case the amplitude of the 
motion is directly limited by the spatial extension of the 
wavepacket coherence. Controlled motion in modulated 
lattices has also been observed experimentally in jH, [l2[ • 
Our work shows how a suitably "overmodulated" driving 
allows to displace a wavepacket along paths of arbitrary 
shape with controlled (and even engineered) dispersion. 

Let us first consider the quantum dynamics of an atom 
in a ID tilted (or "washboard") potential formed by a 
sinusoidal potential superposed to a constant force , or 
a linear potential FxX. We use normalized variables such 
that lengths are measured in units of the lattice step d 



(= Xl/'^^^l = 27r//cL being the laser wavelength) , energy 
in units of the recoil energy Er = h^k\/ (2M) where M 
is the atom mass, and time is measured in units of h/ Er. 
The corresponding Hamiltonian is 



1 

2m* dx^ 



[K + Ax{t)] cos(27rx) + FxX, (1) 



where m* = 7r^/2 is the mass in normalized units, F^ is 
the constant force measured in units of En/d, h = 1, and 
Vx is the lattice amplitude, to which a time-dependent 
component Ax{t) can be added. Tilted optical lattices 
have been experimentally realized by many groups (isl . 

[13, [33. 

Dynamics in a tilted lattice is conveniently described 
by Wannier-Stark (WS) states [H, [3ll-[3i| , which are the 
eigenstates of the Hamiltonian of Eq. ([T]) with Ax{t) = 0. 
In numerical simulations a finite lattice is used, but if 
the bounds of the lattice are far from the region of 
interest there is no essential difference in the dynam- 
ics. The eigenenergies of the Hamiltonian. ([!]) form a 
"ladder" structure separated by the "Bloch frequency" 
ujb = 27r/Te= Fx {ojb = Fxd/h in usual units). De- 
pending on the ratio Vx /Fx each well can host more than 
one WS state, each family of states then forming its own 
ladder j32| . Throughout this paper the parameters of the 
potential and the initial conditions are chosen so that the 
atomic dynamics is accurately described by the lowest- 
ladder WS states only. This situation can be realized 
experimentally by using cold enough atoms and raising 
the optical potential adiabatically pJi]. As WS states 
are, for the parameters used in the present work, highly 
localized on a potential well, we label a WS state by the 
well index n in which it is centered (pn{x). These states 
are invariant under a translation of a integer number n 
of lattice steps, provided the associated energy is also 
shifted by tiub, that is 



En = Eq^ nUJB 



(2) 
(3) 



(in what follows we set £^0 = 0)- It is well known that a 
wavepacket submitted to the Hamiltonian of Eq. ([!]) with 
Ax{t) = has an oscillatory behavior, called Bloch oscil- 
lation, of period T^. Adding a time-dependent potential 
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Figure 1: (Color online) Arbitrary wavepacket motion in two dimensions. The plots display the probability of presence \'^xy{t)f 
integrated on time, and the arrows indicate the sens of the motion, (a) Square path formed of four straight lines, (b) cxo-path, 
obtained by choosing ax = cos (^It), ay = cos (2Qt), /3x — Py — 0. Plot (c) is made is made by combining straight lines and 
arcs. Parameters of the potential are Vx,y = 2.5 and Fx,y = 0.2 (note the different scales in each plot). 



modulated at (or around) the frequency cj^ is thus a good 
way to create a resonant response in the dynamics. 

We can expand the atomic wave function i/jxit) over 
the WS states (of the first ladder): 

V^.W = ^Cn(t)e*^-(^Vn(x), (4) 
n 

where the phase ^n(^) is conveniently defined by: 

t 

(t)n{t) = -nujBt - V^Mo J Ax{t')dt' (5) 



and Mo is the p = value of the coupling parameter M^: 

Mp = {(pi\ cos{27rx) \(pi^p) = {(fo\ cos(27rx) \(pp) . (6) 

[the last identity is a consequence of Eq. (|2])]. Bring- 
ing Eq. dll into the Schrodinger equation associated to 
the Hamiltonian Eq. ([T]) one obtains the following set of 
equations: 

Cn{t) = -lA, (t) M^c^+ne-^^--^ (7) 

Eqs. ([7j) can be simplified by using the fact a WS state 
overlaps significantly only with the WS states associated 
with the nearest neighbor sites, that is Mp ^ {)ioY\p\ > 1. 
We now set in Eq. ([T]) the driving as 

Aa:{t) = {t) sin {cUBt + Pa:) • (8) 

where ax (t) is a slowly varying function (|a~^(iaa;/(it| <C 
ujb)- An additional simplification is obtained by neglect- 
ing fast oscillating terms; one then finds from Eq. ([7]): 

c^t) = ^^^^ [c^-^e-^^^ - w'^^] (9) 



whose general solution is 

Cn{t) = Cn+pit = 0)e"'^^Jp (^-Mi a, {t') dt'^ 

(10) 

{Jpix) is the Bessel function of the first kind and order p) 
which is a straightforward extension of a previous result 
(3H I. The mean position of the wavepacket, averaged over 
a Bloch period T^, is 

(^)t = (^)t=0+ / Vx{t')dt', 

Jo 

where 

Vx{t) = Miax{t)Re{ae-'^-) (11) 

is the instantaneous wavepacket drift velocity and 
a = (0) (0) is the initial coherence of the 

wavepacket. The wavepacket also presents a diffusion 
characterized by an instantaneous diffusion coefficient 
D{t) (X [Miax{t)lm{ae~'^^'')]. Hence, the motion and 
the diffusion can be controlled by changing the temporal 
driving ax and phase /3x • The diffusion can be suppressed 
by setting (3x so that Im (cre~*^^ ) = 0; one then obtains an 
undeformed translation of the wavepacket with a veloc- 
ity given by Eq. ([TT|) . Vx {t) = iMia^^ {t) |<j|. Conversely, 
setting f3x so that Re(cre~*^^) = leads to a purely dif- 
fusive motion with no displacement of the wavepacket 's 
center of mass. 

Consider now a "square" 2D lattice formed by the in- 
terference of two orthogonal pairs of counterpropagating 
laser beams. It can be shown that the resulting Hamil- 
tonian is separable (29|: Hxy = Hx -\- Hy^ where Hx is 
given by Eq. ([Tj) and Hy is the obvious generalization 
for the y coordinate. The solution can be written as 
^^xyit) = ipx{t)'ipy{t) and the 2D dynamics is obtained 
simply by solving two identical ID Schrodinger equations 
for each V^n(^) {u = x, ^), with the corresponding Hamil- 
tonian Hw We can thus induce a controlled 2D motion of 
the wavepacket by choosing suitable lattice modulations 
Au{t) =au{t) sm{ujBt + Pu)- 
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Figure 2: (Color online) The left plot on the top row displays a few wavepackets (or different slices of a large wavepacket) 
arbitrarily distributed in space. The potential defined by Eq. ([13]) sets the different wavepackets in motion with a velocity 
directed towards the origin (situated close to the upright corner) whose amplitude increases with the distance from the origin. 
The slices are thus progressively concentrated around the origin. Plot (b) show the evolution of the total spatial dispersion 
with respect to the origin. Parameters of the potential are Vx,y = 2.5, Fx,y = 0.25 and k/kh — 0.02. 



In order to illustrate the possibilities opened by our 
method, we present in Fig. [1] different trajectories ob- 
tained by numerical integration of the Schrodinger equa- 
tion with fixed parameters and an initial wavepacket of 
gaussian shape 

*xy(0) = 2]exp( )^i{x)^m{y)- (12) 

Note that this form implies that a is real, the "zero diffu- 
sion" condition is thus fulfilled if Px = Py = ^- Plots (a)- 
(c) in Fig. [T] show the square- modulus of the wavepacket 
integrated on time^ and the arrows indicate the sense of 
the motion. Plot (a) shows a square trajectory obtained 
with ax = a^y = 0, for < t < 15Tb , = , = 1 
for 15Te <t< 30Te, and so on. In plot (b), we drive the 
wavepacket into an oo-shaped Lissajous curve by setting 
Oixif) = cosiVtt)^ j3x = ^ and ayit) = cos(21]t), j3y = 0, 
with the over modulation frequency = a;B/250. By 
combining paths one can generate any type of trajecto- 
ries in 2D (and, by an obvious generalization of the above 
discussion, in 3D): Plot (c) in Fig. [T] displays a /3-shaped 
trajectory, where even a turning point has been "drawn". 
In all these numerical experiments the diffusion is very 
small, the width of the wavepacket varies only slightly 
during the evolution. 

New kinds of dynamics can be obtained by using an 
overmodulation whose amplitude varies slowly also in 
space, \.e Ax{t) Ax{x^t) in Eq. ([T]). This can be real- 
ized jl, [13 by adding a second laser beam with a different 
spatial period /c^, which produces a spatial modulation 
of the lattice amplitude corresponding to the beat note 
of the two spatial frequencies. The potential is then 

Au{u^t) = sm{ku) {sin ujBt -\- Pu) (13) 

where k = (/c^ — k'j^)/kL- In the limit \k\ <C 1, siii{ku) ~ 
ku and the velocity of the wavepacket, which is propor- 



tional to the modulation amplitude, varies linearly with 
the position^ as one can deduce from Eq. ([TT]) : Vu = 
Micrku (taking j3u = 0). All parts of the wavepacket, 
whatever their positions, will move towards the origin 
(0, 0), and the closer of the origin it is, the smallest its ve- 
locity. This effect can thus be used to concentrate a large 
wavepacket, or various wavepackets distributed at arbi- 
trary positions, on a few potential wells. Fig. [2](a) shows 
a numerical experiment in which such a concentration is 
realized. Each colored spot represents a wavepacket, or 
a part of a large wavepacket. As times goes by, the dif- 
ferent spots are seen to converge to the same position. 
Plot (b) in Fig. [2] shows the evolution of the total spatial 

1/2 

dispersion ((Ar^) ). This unavoidable dispersion lim- 
its the maximum time during which the overmodulation 
can be applied, and thus the maximum density which can 
be obtained. 

One can also combine temporal and spatial modula- 
tions to produce a rotation of a wavepacket. Consider a 
cigar-shape wavepacket as shown in Fig. [3l^a). Inducing 
a uniform rotation around the origin (0, 0) means to give 
a slice at position r a local velocity v perpendicular to 
r. This can be done by producing a modulation of the 
form 

Ax = - sin{ky) sin {ajst l3x) ^ -ky sin {ajst + /Sx) 
Ay = sin{kx) sin {cjBt -\- f^y) ^ kx sin {uBt -\- /3y) . 

Figure [3] shows an example of wavepacket rotation pro- 
duced by this technique. 

In conclusion, this paper illustrates the almost unlim- 
ited power of overmodulated optical lattices to manip- 
ulate atom wavepackets. The technique is immediately 
applicable to state-of-art experiments, and that can play 
an important role in preparing complex initial state for 
fundamental studies of quantum dynamics, but might 
also useful for playing with atomic qbits in the realm 
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Figure 3: 
text). 



(Color online) Rotation of the wavepacket, obtained by using combined spatial and temporal overmodulations (cf. 



of quantum information experiments. Although we have 
discussed only 2D examples, that are easier to under- 
stand and to display pictorially, all the techniques illus- 
trated above are straightforwardly generalized to the 3D 
case. 
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